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Abstract. We provide a surprising new application of classical approximation 
theory to a fundamental asset-pricing model of mathematical finance. Specif- 
ically, we calculate an analytic value for the correlation coefficient between 
exponential Brownian motion and its time average, and we find the use of di- 
vided differences greatly elucidates formulae, providing a path to several new 
results. As applications, we find that this correlation coefficient is always at 
least l/y/2 and, via the Hermite-Genocchi integral relation, demonstrate that 
all moments of the time average are certain divided differences of the expo- 
nential function. We also prove that these moments agree with the somewhat 
more complex formulae obtained by Oshanin and Yor. 



1. Introduction 
We begin with exponential, or geometric, Brownian motion, defined by 

(1.1) S-lt) = e('-^)*+'^-^(*), t>0, 

where B : [0,c») — > M is Brownian motion. In other words, i? is a stochastic pro- 
cess, or random function, for which B{0) = 1, its increments are independent, and, 
for < s < t, the increment B{t) — B{s) is normaUy distributed with mean zero 
and variance t — s. The basic properties of Brownian motion are explained in Sec- 
tion 37 of Billingsley (1995), while Karatzas and Shreve (1991) is a comprehensive 
treatise. At a more elementary level, Norris (1998) provides a lucid derivation of 
the main properties of Brownian motion, whilst Higham (2004) provides a more 
general introduction well-suited to the numerical analyst. 
We shall study the time average 

(1.2) A{T):=^ S[t)dt 

using the calculus of divided differences, a fundamental tool in approximation the- 
ory. We will, in particular, show that the correlation coefficient between A{T) 
and S{T), the moments of A{T) and, more generally, joint moments of S{T) and 
A{T) can be elegantly, and usefully, expressed in terms of divided differences of 
the exponential function. Now the time average A{T) has been extensively studied 
in the literature of Asian options; see, for instance, Yor (1992), Yor (2001) and 
Oshanin, Mogutov and Moreau (1993). However, we find that our use of divided 
differences both simplify and elucidate formulae. In Section 2, we derive the corre- 
lation coefficient for S{T) and A{T), finding that it is always at least l/-\/2, thus 
explaining the relative high correlation that is observational folklore in the financial 
communityl3 In Section 3, we demonstrate that the divided differences occurring 
in the lower moments of S{T) and A{T) generalise to all moments, using the fact 
that the integral of an exponential function over a simplex can be expressed, via 
the Hermite-Gennocchi formula, as a certain divided difference of the exponential 



"'^I am particularly grateful to my friend and colleague Dr Dirk Siegel, a fellow PhD student of 
Mike Powell, who made this observation known to me when consulting for IBM. 
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function. In Section 4, we provide the divided difference theory required by the 
paper; some of this can, of course, be found in the dedicatee's excellent textbook 
Powell (1981), but derivations of the Hermite-Gennochi and Leibniz formula are 
less easily available. Therefore we have provided their brief derivations in the hope 
that this will enhance the paper's use to both the mathematical finance and the 
numerical analysis communities. Finally, in Section 5, we use our divided difference 
approach to derive a recurrence relation for the moments of A{T). 
We first observe the familiar result 

(1.3) E5(r) = e('-^'/2)Tj£g.TV2z ^ g(r-.V2)Tg<x^T/2 ^ ^rT _ 

Here Z denotes a generic A^(0, 1) Gaussian random variable and we have used the 
standard fact that 



(1.4) Ee^^ = (27r)-i/2 j e^^e-^'/^ dr = (2^)-i/2 J e'^^^^^-^^'-^'Ur ^ e^" '\ 
Similarly, 

E^(r) = T-^ j ES{t)dt 
Jo 



e"^ - 1 

The approximation theorist will immediately recognise the divided difference 
(1.6) EA(r) =exp[0,rT], 

but a sceptical reader might view this as mere coincidence; in fact, it is but the tip of 
an iceberg. We remind the reader that /[ao, oi, . . . , a„] is the highest coefficient of 
the unique polynomial of degree n interpolating / at distinct points oq, . . . , a„ e K, 
which implies /[ap] — f{ao) and 

.r 1 /(oi) - /(Qq) 
fli - ao 

Further, it is evident that a divided difference does not depend on the order in 
which the points oq, ai, . . . , a„ are chosen. As mentioned above. Section 0] collects 
further divided difference theory required by this paper. 

2. The correlation coefficient between the time average and the 

ASSET 

We shall compute the correlation coefficient between S{T) and A{T). Specifi- 
cally, we calculate 

^ E iSiT)A{T)) - E iS{T)) E (AjT)) 
^varS'(r)varv4(r) 
We find an elegant divided difference expression for R. 
Theorem 2.1. The correlation coefficient (|2.ip is given by 

exp[rT, 2rT, (2r + a'^)T] 



(2.2) R = R{rT,a^T) 



^2 exp[2rT, (2r + a'^)T] exp[0, rT, 2rT, (2r + a'^)T] 



Let us begin our derivation. 
Lemma 2.2. IfO<a<b, then 

(2.3) ES{a)S{b)=cxp(^a{r + cr'^) + bry 
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Proof. We have 

ES{a)S{b) = ]E5(a)2p(b-a)(r.-.V2)+-(S(fc)-i3(a)) 

^ g(2r+<T2)Q^(6_a)r 

(2.4) = f,<^{r+<y''))^br^ 

where Z ^ N{0, 1) and we have used (|1.3p . □ 
Proposition 2.3. We have 

(2.5) ES'(T)A(T) = exp[rT, (2r + CT2)r]. 
Proof. Applying Lemma 12.21 we obtain 



ES{T)A{T) = T-^ [ ES{t)S{T) 
Jo 



dt 





= exp[rT,{2r + (t'^)T]. 

Proposition 2.4. 

(2.6) E(A(T)2) =2exp[0,rr, (2r + cr2)r]. 

Proof. We find 



E(A(r)2) = E5(ti)5(t2)dt2) d<i 

(2.7) = 2r-2y" (^J'ES{h)Sit2)dt2)dh. 
Thus 

E(A(T)2) = ^ ' e''(*i+*^)e'^'*Mt2) dii 



2T-2 / e^'^(- ^)dh 

Jo ^ r + cr^ / 

2 

[exp[0, (2r + a^)T] - exp[0, rT]] 



□ 



(r + cr2)T 

(2.8) = 2exp[0,rr, (2r + ct2)T], 

using the divided difference recurrence relation (j4.ip to obtain the final line. □ 

Any reader still doubtful of the simplification provided by divided difference nota- 
tion might consider the alternative expression provided in Hull (2000): 

2g(2r+o-^)T 2/1 

E (A(T)2) = (^_^^2)(2j. + cr2)T2 + Vf2 (,2r + CT2 

There is a similar divided difference relation for K{A{T)™), described in the next 
section, but we now complete our derivation of Theorem 12. II 

Proof of TheoremlKTl Applying (|Lil [T31 1^3)) and (HTTj) . we obtain 

ES{T)A{T) - ES{T)EA{T) = exp[rr, (2r + a^)T] - e''^(e''^ - l)/(rT) 

= exp[rT, (2r + (t'^)T] - exp[rT, 2rT] 

(2.9) = a^Texp[rT,2rT,{2r + a^)T]. 



4 



B. J. C. BAXTER AND R. BRUMMELHUIS 



Further, 
(2.10) 

varS'(r) = E{S{T)^) - {ES{T))^ = e^^''+''">'^ - e^"^^ = cr^Texp[2rT, (2r + cr^)T]. 
and, by (1131 



e — 1 



rT 

= 2 exp[0, rT, (2r + cr2)r] - 2 exp[0, rT, 2rT] 
(2.11) 2cT'^Texp[0,rT,2rT,{2r + a^)T], 

using the divided difference recurrence (j4.1l) once more. Hence 

exp[rT, 2rT, {2r + a'^)T] 



(2.12) iJ: 



a/2 exp[2rr, (2r + (t2)T] exp[0, rT, 2rT, (2r + a^)T] ' 

□ 



It is remarkable that the divided differences appearing in (|2.12l) are coefficients of 
the cubic polynomial interpolating the exponential function at 0, rT, 2rT, {2r+a'^)T. 
We make three further observations: 

(1) Armed with an analytic expression for the correlation coefficient, we can 
apply the exchange option valuation formula of Margrabe (1978) to derive 
the values of certain Asian options, if we are willing to accept that the 
time-average is suitably approximated by exponential Brownian motion. 
We are investigating the numerics of this rather simple approximation at 
present and preliminary results are surprisingly promising. 

(2) The correlation coefficient R{rT, cr^T) is typically close to unity: typical 
values of r , (T and T produce values of R in the 0.8—0.9 range. In fact, we are 
able to prove that the correlation coefficient satisfies R{rT,a^T) > l/-\/2, 
for all r > 0, cr > and T > 0, a surprisingly high lower bound for the 
correlation coefficient. The details of this derivation are too complicated 
to include here, and we refer the reader to Baxter and Fretwell (2007) 
for further details. However, the numerical findings are summarised in 
Figure [H which displays values of the closely related quantity 



exp[a, 2r] exp[a, 2r, r, 0] ' 

for -20 < a < 40 and 0.1 < r < 10. It is easily checked that R{rT, a^T) = 
yJS{rT,{2r + a'^)T)/2, so that the lower bound R>l/\/2 becomes 5" > 1. 
It is plausible that S{r, a) should be a decreasing function of a, for fixed r, 
because correlation should be a decreasing function of volatility. Further, 
it is not difficult to establish the limiting values lima->.-oo S{a, r) — 2 and 
limr_>.oo S{a,r) = 1. However, further analysis is not straightforward, and 
the analysis of Baxter and Fretwell (2007) makes great use of properties of 
divided differences. 

(3) It is natural to ask whether these divided difference expressions are particu- 
lar to exponential Brownian motion. In fact, similar expressions occur when 
exponential Brownian motion is replaced by certain Levy-stable variants: 
see Baxter, Cartea and Fretwell (2007) for further details. 

3. Computing higher moments of A{T) 

We now demonstrate that the neat divided difference formulae obtained for the 
first and second moments of A{T) are not coincidences, but part of a greater pattern 
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Figure 1. S{r, a) for -20 < a < 40 and 0.1 < r < 10 



from which arise new formulae generahsing the moment calculations of Oshanin et 
al. (1993), Yor (1992) and Yor (2001). 
We begin with the iterated integral 

(3.1) EAirr = r dr^ r dTm-1 ■ ■ ■ r dr^ E5(ri) • • • S{t^). 

Jo Jo Jo 

Now, given any point (ri, . . . , t,„) G [0, T]™, let us sort its components into increas- 
ing order, obtaining (ti, . . . ,t„) (say). Then 

ES{Ti)---S{T^)^ESih)---S{trn) 

and 

(3.2) E^(r)'" ^ m\T-"' f dt,n f dt™-i • • • / ' dhES{ti) ■ ■ ■ S{tm). 

Jo Jo Jo 

Our first task is to calculate the integrand, which we complete after a simple lemma. 
Lemma 3.1. For any positive integer k, we have 

(3.3) — 



[S{tf] = exp(fcrt + ^k(k - 1\ 



Proof. This is almost immediate from (jl.4p : 

ES{t)^ — ]Eg'^'('"^'^^)*+'^'=Vt^ — gfe(r-(T^/2)t+<T^fc^t/2 _ gfcrt+CT^tfc(fc-l)/2 

where A^(0,1). 



□ 
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Proposition 3.2. If < ti < t2 < ■ ■ ■ < tm, then 

m 

(3.4) ES{h)S{t2) ■ ■ ■ S{t,n) = exp(^ (r + (m - k)a^) tk 

fc=i 

Proof. Lemma [2.21 comprises the case m = 2. We complete the proof by induction 
on the number of terms m, first observing that, by a standard property of geometric 
Brownian motion, 

(3.5) ES{h)S{h) ■ ■ ■ S(t^) = ES{hr'ES{t2 - ti) • • • 5(t„, - h). 

Applying Lemma l3. II and our induction hypothesis, we obtain 
(3.6) 

m 

'ES{ti)S{t2) ■ ■ ■ S{trn) = exp(mrti+ahim{m-l) (r + (m - e)a'^) {ti - ti)^ 
The ti coefficient in the exponent is given by 



mr 



(m - l)r + a^ti I -m{m - 1) - ^ £ \ = r + a'^ti{m - 1), 



using the elementary fact that m{m — l)/2 = 1 + 2 + • • • + 7n — 1. The coefficients 
of ^2 , ■ • ■ , are as already stated in p.4p . □ 

Thus the desired integral (j3.2p becomes 

E^(T)™ = m!T-" /" dtm f "dt^^i--- [ \tiES{ti)--- S{t^) 
Jo Jo Jo 

(3.7) = m! / dt^ / ••• / dti expiaiti -\ a^t^n), 

Jo Jo Jo 

where 

(3.8) ak = (r + (to — A:)cr^) T, fc = 1, . . . , m. 

The integral displayed in (|3.7p can now be identified as a divided difference using 
a variant form of the Hermite-Genocchi integral relation. 

Theorem 3.3. Let 

(3.9) bk:^kr + a'^k{k~l)/2, fc = 0, 1,.... 
Then 

(3.10) E(A(T))™ =m!exp[6oT,6ir,...,6„T], to > 0. 



Proof Apply Corollarv H31 to ((XT)) and (E^l), using J^l^i k = j{j + l)/2. □ 



The statement of Theorem 13.31 simplifies when r = cr^ , for then the drift term 
in (|l.ip vanishes, that is, we consider S{t) = exp{ay/tB{t)) alone; this is the spe- 
cial case studied by Oshanin et al. (1993) and Yor (1992), for the formulae grow 
much more complicated without the use of divided differences. Therefore we now 
demonstrate that our expression agrees with theirs. 

Theorem 3.4. // we set r = (7^/2 in Theorem \3.!A then we obtain 

E [yl(T)"] = to! exp[0, rT, 2VT, S^rT, . . . , m^rT] 
(3.11) = m\H^^[-m, . . . , -1, 0, 1, . . . , to], 

where Hc{x) :— exp(c^a;^), a; G K, for any positive c. 
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Proof. We simply set r = in Theorem 13.31 and apply (|4.14p . □ 

We can now apply Corollary 14. lOl to derive the formula given in equation (14) of 
Oshanin et al. (1993). 

Theorem 3.5. // we set r = (7^/2, then 

' ^ w^n - (^) (4<-')"C:) +f (T) ^ 

Proof. Applying Corollarv l4.10l to Theorem 13. 4[ we obtain 

'"i'*<^)'"i = ((5s!)i('-^'''"iiCr)<-''''^'"'""''") 

,3.13) - ,,rr' i-li-rr (t) + 1 (T) ■ 

after some straightforward algebraic manipulation. □ 

If we now replace rT by a and m by j in p.l2p . then we obtain equation (14) 
of Oshanin et al. (1993). 



4. Divided difference theory 

Most of the properties of divided differences required here can be found in Chap- 
ter 5 of Powell (1981). However, proofs of the Hermite-Genocchi integral relation 
are less easily available in the Anglophone mathematical literature, as is our par- 
ticular variant of it, although the specialist can find much useful material in the 
treatise of DeVore and Lorentz (1993). We have therefore provided a derivation 
for the convenience of the reader. The Hermite-Genocchi formula and its conse- 
quences are still very much topics of current research; see, for example, Waldron 
(1998). Furthermore, the result is better served in other European languages; see, 
for instance, Gel'fond (1963) for a French translation of a Russian classic, or indeed 
the original Hermite (1878). 

We recall the divided difference recurrence relation. 

Theorem 4.1. 

/[ai,...,a„] -/[ao,...,a„_i] 



(4.1) /[ao,ai, . . . ,a„] 



for any distinct complex numbers ag, . . . , a„. 

Proof. See, for instance, Powell (1981), Theorem 5.3. □ 

If / is sufficiently differentiable, then we can, of course, define divided differences 
for coincident points. Further, the elementary relation 

(4.2) /[ao,ai] ^ ~ ^ / f ((1 ~ t)aQ + tai) dt, whenacOiGM, 

can be generalised to obtain the Hermite-Genocchi formula. 
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Theorem 4.2 (Hermite-Genocchi) . Let f e C(")(M) and let ao, ai, . . . , a„ be (not 
necessarily distinct) real numbers Then, for n > 1, 

f[ao,ai, . . . ,a„] 

/*^"^(toao + hai H h t„a„) dti-- ■ dtn, 



's„ 

(4.3) = dti dt2--- dtnP'''{toao + tiai + ---+tnan) 

Jo JO Jo 

where the domain of integration is the simplex 

(4.4) Sn = l^t={ti,t2,...,tn)eWl: J2tk < l| 
and 

n 

to = I — ^ tk- 

k=l 

Proof. We shall prove (|4.3p by induction on n, observing that 

f'{toao + tiai)dti^ [ /'(oq + <i(ai - aq)) dti = /(«o) ^ f[ao^ai] 

Si Jo 0.1 ~ «o 

To extend the formula to higher-order divided differences, we note that 
(4.5) 

„r /[ai,a2,...,a„+i] -/[ao,ai,...,a„] 
j[ao,ai, . . . ,a„,a„+ij — — g[ao, On+i], 

In+l — 0,0 

where 

(4.6) g{x) = f[ai,...,an,x], x&R. 
Now 

g{x) = [ {xto + aiti H h a„i„) dti ■ ■ ■ dtn 

so that 

g'ix)^ I tof^-+'\xto + fliti + • • • + a„t„) dii • • • dtr, 



Therefore 

/[ao, Oi, . . . , a„, On+i] 



1 

dr .g'((l - T)ao + Ta„+i) 



dr / dti--- dtn io/^"+^H[(l - T)ao + Ta„+i]to + ai<i + 
/ dti--- dtn / dr io/^"+'H[(l - T)ioao + V a^te + rtoOn+i) 

JS„ JO 

1 /-i-ti /-i-Efcii'tfe / »+i 



dii di2 • • • di„+i/("+i) ( roao + t^a,. 1 

j(n+i) (j-gap + tiai H h tn+ian+i) dti - - - di„+i, 

-1 

where we have used the substitution tn+i — tor and the notation To — 1 — J2k=i ^k- 

□ 
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We shall need a variant form of the Hermite-Genocchi integral relation for which 
the following notation is useful. Given any real n x n nonsingular matrix V, with 
columns wi, . . . , w„, we let K(V) denote the closed convex hull of 0, wi, . . . , u„, i.e. 

K{V) := conv{0, vi, . . . , u„}. 

In this notation, the Hermite-Gnocchi integral relation states that 

(4.7) /[ao,ai, . . . ,a„] = / /^"M ao + (a - aoe)'^ y ) dy, 

where 

















) 




Van/ 







and /„ denotes the n x n identity matrix. Integrating the nth derivative over the 
simplex K{V) yields a useful variant form of Hermite-Genocchi. 

Theorem 4.3. Let V E M"^" be any nonsingular matrix. Then 
(4.8) / (a^y) dy = /[O, (V^a),, . . . , (V^aU 

where {V'^a)^ denotes the kth component of the vector V'^ a. 
Proof. Substituting y — Vz, Hermite-Genocchi implies the relation 

{(V^afz) dz = /[O, (y^a)i, . . . , (F^a)„]. 



K(V) 



□ 



Corollary 4.4. For any function f G C^"^(]R), we have 

/ dXn I dXn-1--- / dxif^"^ I y^flfcXfc ) 

Jo Jo Jo Vfc=i / 
(4.9) = /[0,a„, a„ + a„_i, . . . ,a„ + a„_i H h ai]. 

Proof Set 

V 



/I 

1 1 



\1 1 ... Ij 

in Theorem 1131 □ 



The exponential function is a particularly important case for us, in which case 
the Hermite-Genocchi formula becomes 

(4.10) exp[ao,ai,...,a„] = / e*''°''+*i"i+-+*"°" dti . . • dt„ 

and Corollarv 14.41 takes the following form. 
Corollary 4.5. We have 

/ dXn / dXn-1 ■ ■ ■ I dxi exp ^ flfeXfe 

Jo Jo Jo Vfc^i y 

(4.11) = exp[0,a„, a„ + a„_i, . . . ,a„ + a„_i H h ai]. 
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Proof. Let / be the exponential function in Corollary 14.41 □ 

Further, we note that, for the exponential function, Theorem 14.31 becomes the 
interesting formula 

(4.12) tAtTi f e-"ydy = exp[0,{V^a)i,...,{V^aU. 
|aet V I Jk(v) 

Thus, integrating exponentials over simplexes or, more generally, a polyhedron 
formed by the disjoint union of simplexes, will generate divided differences of the 
exponential. 

We shall also need two simple preliminary results. Let us use P„ to denote the 
vector space of polynomials of degree n. 

Lemma 4.6. We have 

(4.13) exp(^) exp[Ao, . . . , A™] = exp[Ao + /i, . . . , A™ + /x], 
where Aq, . . . , A,„ and /i can be any complex numbers. 

Proof. Immediate. □ 

Lemma 4.7. Let / : C — >■ C and /e< ai, . . . , a„ be distinct nonzero complex numbers. 
Then 

(4.14) /[O, aj,..., al] = g[-an, -oi, 0, ai, . . . , a„], 
where g{z) ~ /(^^); for z € C. 

Proof. Let p G P„ interpolate / at 0, af, . . . , a^. Then q{z) := piz"^) is a polynomial 
of degree 2n satisfying q{±aj) = p{a^j) = f{aj) = g{±aj), for j = 0, . . . , n, setting 
ao = 0, for convenience. The result then follows from uniqueness of the interpolating 
polynomial. 

□ 

It is well-known that a divided difference at equally spaced points can be ex- 
pressed in a particularly simple form using the forward difference operator 

A,/(x) ■.= f{x + h)-f{x), 

which wc shall need when demonstrating the equivalence between our moment 
calculations and those of Oshanin et al. (1993) and Yor (1992). The next proposition 
is well-known, but we again include its short proof for the reader's convenience. 

Proposition 4.8. Let / : M — > M, let h be any positive constant and let n be a 
non-negative integer. Then 



(4.15) f[x,x + h,x + 2h, . . . ,x + nh]^ 



i\h" 



Proof. It is easily checked that f[x,x + h] = Afif{x)/h. Further, if we assume 
(|4.15p for n — 1 , then the divided difference recurrence relation implies that 

f[x, X + h, . . . ,x + nh] 

f[x + h, . . . , X + nh] — f[x, X + h, . . . , X + (n — l)h] 
nh 

^ Ahf[x,...,x + {n-l)h] 
nh 

1 Arv(^) 



nh \{n-iy.h''-^ 



FUNCTIONALS OF EXPONENTIAL BROWNIAN MOTION AND DIVIDED DIFFERENCESll 



Thus the rcsuh fohows by induction. □ 
Corollary 4.9. Let / : K ^ M and let h be any positive constant. Then 

(4.16) fix, x + h,x + 2h,...,x + nh] = -—Y,[l) (-l)""'/(a: + kh). 



Proof. We define the forward shift operator 

Ehfix) f{x + h), xe R, 
and observe that, by the binomial theorem, 

Alfix) = - ir fix) - (l) {-ir-^Eif{x) = it] {-ir-'f{x+kh). 



k=0 ^ ^ fc=0 ^ 



□ 



Corollary 4.10. Let / : R — > M and let h be any positive number. Then 
(4.17) 

f[-nh, -{n - l)h, ...,-h,0,h,..., nh] = — -— EL" (^^^'^^^ " 

Proof. This is an immediate consequence of Coroharv 14.91 □ 

We shall also need the Leibniz relation for divided differences of a product when 
deriving the recurrence differential equation for moments. 

Theorem 4.11 (Leibniz). Let D be any subset of C containing the distinct points 
zq, zi, . . . , Zn and let v and w be complex-valued functions on D. Lf u = v ■ w, then 

n 

(4.18) u[zo,...,Zn] = y^^v[zQ, . . . , Zk]w[zk, . . . , Zn]. 

k=0 

Proof. Let p G P„ be the unique polynomial interpolant for u written in standard 
Newton form, that is, 

(4.19) p{z) = v[zo] + v[zo,zi]{z - zo) H v[zo,zi, . . . ,z„](z - zo) • • • (z - z„_i). 

We shall let g G P„ be the unique polynomial interpolating w, but with the points 
chosen in the order z„, z„_i, . . . , zq, that is, 

(4.20) q{z) = w[zn] + w[zn, z„_i](z - Zn) 4 h w[z„, . . . , zo](z - z„) • • • (z - Zi). 

Now their product p • g is a polynomial of degree 2n. Dividing this polynomial by 
(z — zo) • • • (z — Zn), we obtain 

p{z)q{z) = r(z) + s(z)(z - zq) • • • (z - z„), 

where r G Pn- We see that u{zj) — v{zj)w{zj) — p{zj)q(zj) — r{zj), for < j < n. 
Hence, by uniqueness of the polynomial interpolant for u in P„, we obtain 

(4.21) r(z) = u[zo] H h u[zq, . . . , z„](z - zq) • • • (z - z„). 

We obtain (|4.18p by equating the coefficients of z" in (14.2ip and the product of the 
expressions in (|4.19|) and (|4.20p . modulo (z — zg) • • • (z — z„). □ 
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5. A RECURRENCE RELATION 

The Feynman-Kac formula (see, for example, Karatzas and Shreve (1991)) sug- 
gests that the moments -E„(t) := E,{A{t)"-) of the time average should satisfy a 
certain differential equation, which we shall also obtain as an illustration of the 
divided difference approach. 

Theorem 5.1. Let {c„}J^]^ be any strictly increasing sequence of positive numbers 
and define e„ : (0, oo) — >■ R 6?/ the divided difference 

(5.1) e„(i) = exp[0, cit, . . . , c„i], i > 0, n>0. 
Then 

(5.2) te'Jt)^enit){cnt-n) + en-lit), forn>l. 



Proof. Applying the Hermite-Genocchi formula, we obtain 
(5.3) e„{t) = I exp[tc^y) dy, 



where b — (ci, . . . , Cn) , and differentiating (j5.3p yields 

(5.4) e'n{t)= [ eMtc^y){c^v)dy. 

Now writing g{s) = s and applying Leibniz's formula for divided differences, we 
find 

{g ■ exp) [0, cit, . . . , Cnt] = g[0] exp[0, cit, . . . , c„t] + ^[0, cit] exp[cit, . . . , c„t] 

(5.5) = exp[cit, . . . , c„i]. 

Further, the relation [g ■ exp)*-") = g ■ exp +n exp and (|5.5|) imply 

te'nit) = / (5 • exp)^"^ [tc^v) dy~n exp{tc^y) dy 

= [g ■ exp) [0, cii, . . . , c„t] - n exp[0, cit, . . . , c„t] 

(5.6) = exp[cit, . . . ,c„t] - ne„(i). 

However, 

exp[cii, . . . , Cnt] ~ en-i{t) 

(5.7) en{t) = , 

Cril 

by the divided difference recurrence relation, so that 

(5.8) exp[cit, . . . , Cnt] = CntCnit) + e„_i(t). 

Substituting in provides ([O]). □ 

The corresponding differential equation for En is now immediate. 
Corollary 5.2. The moments satisfy 

(5.9) tE'^{t)=En{t){bnt-n)+En-i{t), forn>l, 
where 6„ is given by p.9p . 

Proof. We apply Theorem 15.11 and Theorem 13.31 □ 
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